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ABSTRACT 

We combine the high-resolution Aquarius simulations with three-dimensional models of reionization 
based on the initial density field of the Aquarius parent simulation, Millennium-II, to study the impact 
of patchy reionization on the faint satellite population of Milky Way halos. Because the Aquarius 
suite consists of zoom-in simulations of halos in the Millennium-II volume, we follow the formation 
of substructure and the growth of reionization bubbles due to the larger environment simultaneously, 
and thereby determine the reionization redshifts of satellite candidates. We do this for four different 
reionization models, and also compare results to instantaneous reionization. Using a simple procedure 
for selecting satellites and assigning luminosities in the simulations, we compare the resulting satellite 
populations. We find that the overall number of satellites depends sensitively on the reionization 
model, with a factor of 3-4 variation between the four models for a given host halo, although the 
difference is entirely in the population of faint satellites (My > —10). In addition, we find that for 
a given reionization model the total number of satellites differs by 10-20% between the patchy and 
homogeneous scenarios, provided that the redshift is chosen appropriately for the instantaneous case. 
However, the halo-halo scatter from the six Aquarius halos is large, up to a factor of 2-3, and so is 
comparable to the difference between reionization scenarios. In order to use the population of faint 
dwarf galaxies around the Milky Way as a probe of the local reionization history, then, it is necessary 
to first better understand the general distribution of substructure around Milky Way-mass halos. 
Subject headings: cosmology: theory - galaxies: dwarf - early universe - methods: numerical 
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1. INTRODUCTION 

A major challenge to the ACDM cosmological model is 
its over-prediction of substructure on small scales com- 
pared to the actual observed number of dwarf galax- 
ies in the Loc al Group (th e so-c alled "missing satel- 
lite problem"; iMoore et all 119991 : iKlvpin et aiTll999f ). 
The most well-studied explanation for this discrepancy 
is that the lowest mass halos were inefficient at form- 
ing stars, so that most of these substructures remain 
dark. Indeed, several processes can suppress star for- 
mation and reduce the number of vi sible halos, in- 
cludi ng tidal stripping of satellites (e.g.. iKravtsov et"all 
[200l . supernova feedback, and the increase of gas tem- 
peratur e due t o cosmic reionization (lEfs tathiou 1992; 
GnedinI 120 00": Be nson et all l2002al: ISMnervillel 120021: 
Madau et al . 2008.) . Dynam ical interactions between 



dwarfs ( D'Onghia et al.ll2009| ) or between dwarfs an d the 
luminous disks of large halos (|D'Onghia et al.l[2010( ) may 
also play a role. 

The discovery of a population of "ultra-faint" dwarf 
galaxies (UFDs; L < lO^Lp,) in the Sloan Dig- 
ital Sky Survey (SPSS) Ce.g.. I\yillman et all [20051 
iBelokurov et al.l 120061 : iZucker et al.l l2006( ) has roughlv 
doubled the number of known satellites, alleviating the 
problem somewhat. Based on the sky coverage and sen- 
sitivity of SDSS, estimates of a complet eness-corrected 
luminosity function down to My ~ —2 (|Koposov et al.l 
120081 iToilerud et all I2008D are available. The possible 
number of faint satellites could be up to the order of hun- 
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dreds, but due to the incomplete sky coverage of SDSS, 
the exact number remains uncertain. 

The discovery of these low-luminosity dwarfs has also 
raised many new questions, however. Follow-up ob- 
servations have shown that they are the faintest, most 
dark matter-dominat e d and metal-poor systems known 
(jSimon fc Gehall2007t IKirbv et al.l[200l . with metallic- 
ity pattern s similar to the low est metallicity stars in the 
halo (e.g.. IFrebel et al.ll2010D . In order to understand 
the cosmological context of UFDs, including their rela- 
tionship to the first galaxies, and to the building blocks 
of the stellar halo, it is important to be able to model 
condi tions in the universe at the time of their formation 
(e.g., IFrebel fc Bro^[20Tol ). 

Semi-analytic models of galaxy formation have been 
able to explain the observed number of satellites by 
invoking various feedback mech anisms, inc luding an 
external reionizat i on field (e.g. , iBenson eF al. 20021^ 
IMaccio et aIll2010l:ILi et al.ll2010( ). According to this sce- 
nario, as the intergalactic medium is heated by reioniza- 
tion, gas can no longer be accreted by and cool within 
halos below a certain virial temperature, inhibiting fur- 
ther star formation in the lowest mass halos. How- 
ever, these models generally rely on a simplified, uni- 
form reionization field, which is likely quite different 
from reality. Simulations of reionization indicate that 
it was likely extended in time and spatially inhomoge- 
neous, with a wide range in re ionization redshifts for 
different parts o f the universe (jBarkana & Loeb '200i 
, Lidz et al.. ,20071: iMcQuinn et al.li2007; .Ahn et al.,,2009 : 
Ulvarez et al.ll2009ir 
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Several recent studies have examined various aspects 
of inhomogeneous (patchy) r eionization and its conse- 
quences for galaxy formation. iMunoz et alj ()2009D com- 
bined a reionization model with merger trees drawn 
from the Via Lactea II simulation, and used the ob- 
served pro£erties_of satelli tes to constrain model param- 
eters. lAlvarez et al.l ()2009[ ) combined N-body and three- 
dimensional reionization calculations to investigate the 
relationship between rei o nizat ion history and local en- 
vironment. iBusha et al.l (|201 (/) used their resulting dis- 
tribution of reionization redshifts together with the Via 
Lactea II subhalos, and found that the resulting number 
of Milky Way satellites can vary by an order of mag- 
nitude. Understanding the connection between environ- 
ment, reionization epoch, and faint satellites, then, can 
in principle help us to understand how the Milky Way 
formed. 

Our work is similar in spirit to the above efforts, in 
that we aim to examine how the history of reionization 
impacts the abundance of faint satellites around Milky 
Way-like halos. However, our methods differ from the 
earlier studies in that we combine reionization calcula- 
tions with high resolution zoom-in simulations of Milky 
Way-sized halos in the same box, enabling us to track the 
growth of reionization bubbles and dark matter substruc- 
ture simultaneously. This means that, subject to the 
approximations inherent to our treatment, the reioniza- 
tion history is completely determined by the large-scale 
environment, and infalling subhalos can have different 
reionization redshifts than the Milky Way-like final host. 
Moreover, since we have six high resolution halos avail- 
able, we can roughly estimate the consequences of cosmic 
variance, rather than assuming that a single dark mat- 
ter halo is representative of the Milky Way. Comparing 
the results from the patchy scenarios with instantaneous 
and homogeneous reionization models, we aim to quan- 
tify how much of a difference ignoring the patchiness of 
reionization makes, and to what extent the faint satellites 
trace the local reionization epoch. 

This paper is organized as follows: Section [2] describes 
the simulations that are the basis of our analysis: the 
Aquarius suite of Milky Way halos, and reionization cal- 
culations based on its "parent" simulation. Millennium 
II. Our approach is described in Section [31 detailing how 
we pick out satellite candidates from the merger trees, 
and assign luminosities to them. Section |4] describes our 
results, showing luminosity functions in different scenar- 
ios, and exploring the differences between the patchy and 
homogeneous models. Our conclusions are summarized 
in Section [5j 

2. SIMULATION DATA 

2.1. The Aquarius simulations 

Our modeling is b ased on the AquariufQ simulations 
(jSpringel et al.ll2"008[ ) , a suite of six highly resolved Milky 
Way sized dark matter halos. These are re-simulations 
of six halos from a lower resolution version (900'^ par- 
ticles) of the Millennium-II simulation, which is a peri- 
odic box of size 100/i~ ^ Mpc containing 2160^ particles 
(|Bovlan-Kolchin et alSToOQh . The six halos were selected 
on the basis of final mass, and by not having a massive 

^ http:/ /www. mpa-garching.mpg.de/aquarius/ 



Table 1 

Basic parameters of the six Aquarius halos. 



Name 


mp (Mq) 


M200 (Mq) 


r200 (kpc) 


Aq-A-2 


1.370 X 10"' 


1.842 X IQi^ 


245.88 


Aq-B-2 


6.447 X 10^ 


8.194 X 10" 


187.70 


Aq-C-2 


1.399 X 10'' 


1.774 X IOI2 


242.82 


Aq-D-2 


1.397 X 10'' 


1.774 X I012 


242.85 


Aq-E-2 


9.593 X 103 


1.185 X I012 


212.28 


Aq-F-2 


6.776 X 103 


1.135 X I012 


209.21 



close neighbor at z = (no late-time major merger); 
otherwise they were selected randomly. We note that 
one of these halos (the F halo) is somewhat of an outlier 
in terms of its merger history, in that it experienced its 
last major merger at z = 0.7. 

The same cosmological parameters wer e adopted as for 
the original Millennium simulation ( Sprin ge! et al.i r2005'): 
n,n = 0.25, r^A = 0.75, as ^ 0.9, "n^ = 1 and Hub- 
ble constant Hq = 100ft- km s^^ Mpc^^ — 73 km s~^ 
Mpc"^ While these are consistent within the uncertain- 
ties with the parameters estimated from the three-year 
Wilkinson Microwav e Anisotropy Probe (WMAP) data 
(jSpergel et al.l 120071 ). the value of cts = 0.9 is slightly 
higher than that inferred fro m the most recent measure- 
ments (jKomatsu et al.l [201 It ). In principle, the adopted 
as could lead to an overestimate of the number of halos 
collapsing around the reionization epoch, but studies like 
iBovlan-Kolchin et al.l (|2011aD suggest this would have a 
relatively minor effect. 

All six halos (labeled 'A' to 'F') were simulated with 
at least two different resolutions; we use the higher of the 
two ("level 2"; softening length of 65 pc) in our analysis. 
Table [T] lists the particle mass, final halo mass, and virial 
radius of each of the halos. Mass and radius are listed as 
Af2oo and r2oo, defined as the mass enclosed in a sphere 
with mean density 200 times the critical density, and the 
corresponding virial radius. 

Snapshots are stored at 128 output times, from red- 
shift 127 to 0, equally spaced in log a where a = 1/(1 -I- z). 
For each snapshot, subhalos are identified using the SUB- 
FIND algorithm (Springcl et al. 2001), and linked be- 
tween snapshots in a merger tree. Since SUBFIND re- 
quires a minimum of 20 gravitationally bound particles 
to define a subhalo, the smallest halo masses we resolve 
are of order a few times 10^ solar masses. The merger 
trees are our starting point for identifying subhalos in the 
simulation that may host present-day luminous satellite 
galaxies. 

2.2. Reionization calculations 

The Aquarius halos are drawn from a larger parent sim- 
ulation, so we can use the parent box to calculate how 
regions around these halos reionize depending on their 
en vironment. We u se the semi-analytic approximation 
of lZahn e^ aT (200 71), based on t he excursion set treat- 
ment of iFurlanetto et al.l (|2004bl [al). This semi-analytic 
model assumes that galaxies reside in halos above some 
minimum mass, Mmin, and that the number of ionizing 
photons produced by a galaxy is directly proportional to 
its host halo mass. A given region is considered ionized 
when the collapse fraction - i.e., the fraction of mat- 
ter that lies in halos above the minimum mass - in the 
region exceeds some threshold value. This threshold de- 
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Table 2 

Properties of the reionization models 



Model No. 


^reion ^ 




^end 


Median 


Te 


1 


9.31 


1.44 


7.44 


8.84 


0.075 


2 


10.71 


1.35 


8.92 


10.27 


0.091 


3 


11.97 


1.28 


10.27 


11.57 


0.108 


4 


13.50 


1.20 


11.88 


13.14 


0.129 




Figure 1. Slice through a reionization map from Model 1, illus- 
trating the patchy nature of the process. The color indicates the 
redshift of reionization. The slice is 100 Mpc on a side, and 
195 fe-l kpc (the cell size in the reionization calculations) thick. 

pends on the efficiency of the ionizing sources, and is 
smaller for models with more efficient sources. For the 
purpose of modeling reionization, we fix the minimum 
mass at Mmin = lO^M©, and consider several models for 
the source efficiencies, spanning a range of possibilities 
for the timing and duration of reionization. 

We take the initial conditions of the Millennium-II sim- 
ulation as input to generate reionization maps for the 
cube of 100 Mpc per side. Reionization redshifts are 
tagged spatially, at a resolution of 512 cells per side, or 
about 195 kpc. As an example, a slice through one 
of the reionization cubes is shown in Fig. [1] illustrating 
the patchiness of the processQ 

We calculate four different models, including a very 
efficient reionization that completes around redshift 12, 
and three gradually more extended models. The distribu- 
tions of reionization redshifts in each model, over the cells 
in the box, are summarized in Tableland Fig.[2j In the 
Table, < Zreion > and CTz refer to the mean and standard 
deviation of the reionization redshifts across the cells in 
the box, and z^nd is the redshift where the process is com- 
plete. The quantity is the Thomson optical depth; 
the 7-year WMAP va lue for this is = 0.088 ± 0.015 
(jKomatsu et al.ll20lH ). 

With these maps, we can then find the reionization red- 
shift of any given subhalo as follows. At each snapshot, 
we locate the position of the subhalo in the parent box, 
and look up the redshift at which this region of the par- 
ent box reionizes. The first snapshot where this redshift 
is earlier than the redshift of the snapshot, flags the sub- 
halo as having entered a reionized region, and we define 
this as the reionization redshift of the subhalo. Note that 
this allows not only for different reionization redshifts for 
the six main Aquarius halos, but also that their subhalos 
ending up as satellites are not required to reionize at the 
same time as the host. 

3. METHODOLOGY 
3.1. Picking out satellite candidates 

^ The reionization maps are available from the first author upon 
request. 




Figure 2. Boxplot summarizing the distributions of reionization 
redshifts over the cells in the 100 Mpc box, for the four dif- 
ferent models we considered. The boxes enclose the interquartile 
range, and the horizontal line is the median in each case. The 
whiskers extend to the highest and lowest value still within 1.5 X 
the interquartile range, and the 2nd, 9th, 91st and 98th percentile 
are marked by diamonds. 

We use the level-2 merger trees of the six Aquarius 
halos to determine which dark matter halos can host lu- 
minous satellite galaxies based on several criteria. These 
again reflect our assumptions about when a subhalo is 
able to accrete gas from the IGM, cool its gas, an d form 
stars, f ollowing methods deve l oped in e.g. Madau 'et alJ 
(pOOSl ): iKoposov et all (pOOl ): iBusha et al.l (|20101. We 
only consider subhalos that end up within the virial ra- 
dius of their host halos (see Table [T]). 

First, we require that subhalos must grow larger than 
a minimum Vmax of 15 km/s before star formation can 
take place, and also that this happens before the sub- 
halo is reionized. This corresponds roughly to the virial 
temperature where atomic hydrogen is no longer able to 
cool lO"* K). We do not explicitly include H2-cooling 
"minihalos" as sites for low-mass star formation before 
reionization, but recognize that such halos may have 
been the sites of the very first star formation, whose 
metals see ded the next generation of lower-mass, long- 
lived stars (iBromm et al..,1999. 200llYoshida et alJl2004l : 



iBromm et al.ll2009D . 



Once a subhalo enters a reionized region (as described 
in Scction r2.2l) . the minimum Vmax needed to sustain star 
formation increases, reflecting both the increased tem- 
perature of the IGM to a few times lO'* K due to pho- 
toionization, and the possibility of photoeva.poration of 
gas g ut of smal l halos (jBarkana fc Loeblll999tlliiev et al.l 
[2005h . IGnedinI ([200l found the filtering mass to corre- 
spond to the scale at which the gas fraction in halos 
is significantly reduced due to reionization. Gnedin's 
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expression, however, does not take into account that 
only the cold gas r emai ns available for star formation; 
iMuiioz et all ()2009[ ) and iBusha et all ()2010D argue that 
only halos that grow to a viral temperature ~ 10^ K 
will be able to sustain star formation after reioniza- 
tion. Accordingly, we only let star formation continue 
post-reionization in subhalos whose Vmax is greater than 
50 km . This represents a rather abrupt reionization 
effect; we also run a model where this suppression is 
weaker, with the threshold set at 30 km s^^. 

In addition, we assume that once a satellite starts in- 
teracting with the main halo, its gas is stripped and no 
further star formation takes place. Thus, there are two 
important redshifts associated with each satellite: z^e 
and Ziiifaii, the times when it is reionized and falls into 
the main halo, respectively. 

In order to account for tidal stripping of satel- 
lites, which has been proposed as another important 
mechanism to resolv e the missing satellite problem 
(jKravtsov et al.|[2004[) . we tag the 1% most bound par- 
ticles of the subhalos as star-carrying when they are ac- 
creted onto the main halo. We then track all these tagged 
particles to redshift zero, and determine the final lumi- 
nosity of a satellite based on its remaining tagged parti- 
cles. In this way, we are able to take into account tidal 
stripping of satellites and other dynamical effects, under 
the assumption that the stars are found in the bottom of 
the potential well. 

Of course, a number of infalling systems will be com- 
pletely disrupted, and contribute to the accreted stellar 
halo. In Fig. [31 we show the projected distribution of 
all the dark matter particles that were tagged as stars 
on infall, in reionization scenario 1. Surviving satellites 
here appear as bright, concentrated dots, while the parti- 
cles from the shredded systems make up the diffuse halo. 
Since the accreted stellar halos and streams forming in 
th e Aquarius s imulations have been studied extensively 
in iCooper eTa l. (2010) and Hclmi ct al. (2011), we wiU 
not discuss them further here. We note, however, that al- 
though our model is very simple, it still reproduces mo st 
of the features seen in Fig. 6 of ICooner etlll (1201(11 ) - 
we take this as justification that our simple model still 
reasonably captures the relevant processes. 

We determine satellite candidates for all four of our 
reionization calculations. In addition, we do the same 
for four "instantaneous" redshift models, where we set 
the reionization redshift e qual to the me a n reio nization 
redshift in the entire box. I Alvarez et al.l ()2009[ ) showed 
that the distribution of reionization redshifts for Milky 
Way-sized halos is nearly similar to the global spatial dis- 
tribution, making the global mean an appropriate choice. 
Comparing the resulting satellite populations from the 
patchy and instantaneous reionization models then al- 
lows us to isolate any effects of reionization not being 
spatially uniform. 

3.2. Assigning luminosities to halos 

In order to compare our satellite candidates to the ob- 
served Milky Way population, we need to assign lumi- 
nosities to the subhalos. Ther e are many different pro- 
cedures for doing this; see e.g., iKoposov et all ([20091 ) for 
a discussion. Here, we adopt a simple model where the 
star formation rate is proportional to the mass, that is 



Mstar = / £ X fgas X Mj) M dt, whcrc c = 1 / T IS & constant 
that sets the star formation efficiency (i.e., r can be con- 
sidered a timescale for which a subhalo would convert all 
of its gas into stars). We take /gas xe = 0.08 x 10~^°yr~^, 
tuned both so that we get a reasonable satellite popula- 
tion, but also so that the mass of the accreted stellar halo 
(i.e., the total mass of satellites that fall in and are com- 
pletely shre dded) falls in a mass range around ^ IO^Mq 
(jBell et al.l[2008l : lTumlin son 2010). This second criterion 
is insensitive to the reionization prescription, because the 
mass of the stellar halos is dominated by a few signifi- 
cant progenitors, while the smaller systems affected by 
reionization contribute a small fraction of the total mass. 
While in principle we could increase this product to get 
a better fit for the more extreme reionization scenarios, 
this would then lead to very massive stellar halos. 

Our model is deliberately simple because we wish to 
assess and compare the impact of different reionization 
assumptions without getting addressing uncertainties in 
the baryonic physics. While we do not explicitly model 
effects like the impact of supernova feedback on the satel- 
lites, we assume that the typical effect is captured by set- 
ting e (i.e., the star formation efficiency) appropriately. 
(For an example of a more extensive approach in mod- 
eling the baryonic content of dark halos in cosmolo gical 
simulations, we refer the reader to lFont et al.ll20lTl ) 

Having determined the stellar mass formed, we then 
assign each stellar population an associated present- 
day luminosity, using the stellar po p ulatio n synthesis 
(SPS) models of IBruzual fc CharlotI (12001 . assuming 
a Salpeter initial mass function, and a metallicity of 
Z = 0.0001 (the lowest considered in their model, and 
which captures the typical range of t he ultra-faint satel- 
lites we are most intereste d in, e.g., iKlrbv et all 120081 : 
iSimon et al] 120101: iNorris et~all I2O1O0 . We then track 
the particles tagged with stellar populations to redshift 
zero, and so determine the present-day luminosities of 
the satellite candidates picked out as described earlier. 

As an alternative to the simple parametric method de- 
scribed above, we al so try the abundanc e matching tech- 
nique developed by IBusha et al.l (|2010[ ) to assign lumi- 
nosities t o dark matter halos. T his is an extrapolation of 
the fit of iBlanton et al.l (|2005j ) based on galaxies in the 
SDSS (down to = —12.375) to fainter magnitudes, 
by a power law: 



51og/i = 18.2 - 2.5 log 



1km s 



7.1 



(1) 



The appeal of the abundance matching procedure is that 
there are no parameters to tune, using a model derived 
from statistics of observed galaxies. While i t is an extrap- 
olation at UFD luminosities, .Busha et al.l (j2010t ) showed 
that when applied to the Via Lactea subhalos, the abun- 
dance matching model did well i n reproducing the ob- 
served luminosity function (but see lBovlan-Kolchin et al.l 
I2011bl for a discussion of potential problems with abun- 
dance matching for the MW dwarfs). For purposes of 
comparison, then, we also separately assign magnitudes 
to our subhalos using eqn. [TJ In order to select an ap- 
propriate Vmax for the calculation, we use the peak Vmax 
of the subhalos during the time they were able to form 
stars, as defined in Section l3Tl For most subhalos, this 
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Figure 3. Projected distribution of all dark matter particles that are tagged with stellar populations from infalling subhalos, at z = 0. 
The scale is in kiloparsecs. Both the concentrated, bright satellites and the shredded satellites that make up the diffuse accreted stellar 
halo and streams are visible. 



either corresponds to the time of reionization or time of 
accretion. 

As some basic tests of our model, we check how weU it 
reproduces various features of the observed Milky Way 
dwarf population, such as their radial distribution, and 
observed mass-to-light ratios. Fig. 2] shows the radial 
distribution of satellites in our simulation, as a fraction 
of total number. Wc find all six halos to have a satellite 
distribution consistent with what is observed around the 
Milky Way, with the possible exception of the B halo, 
which has a larger fraction of satellites somewhat closer 
in. 

As for mass-to-light ratios, we note that subhalo 
masses in a simulation cannot be directly compared to 
observed dw arf g alaxy mass e s. H owever, iWalker et al.l 
((20091 IMS) and IWolf et al.l (I2OIOI) show that the mass 
within the half-light radius is well constrained for the 
Milky Way dwarfs. We thus estimate the half-light ra- 
dius of our satellites from their assigned lu minosities by 
fitting a power law to the data tabulated in Walke r et al.l 
((Ml . and measure the m ass within the half-lig ht radius 



by fitting an NFW profile ((Navarro et al.lllQQTl ). The re- 
sult is shown in Fig.[S] the red points are simulated satel- 
lites, and the black crosses (UFDs) and triangles (clas- 
sical dSphs) show the data from the Milky Way dwarf 




200 
r (kpc) 



400 



Figure 4. Radial distribution of satellites, as a fraction of total 
number. The points with error bars show the observed Milky Way 
distribution, where non-SDSS satellites are weighted by a factor of 
fDR5 = 0.194 of the total to account for incomplete sky coverage. 
The theoretical estimates are for reionization Model 1 with strong 
suppression, but the trends are similar for the other models. In 
particular, all halos have a radial satellite distribution consistent 
with that observed around the Milky Way, with the possible ex- 
ception of the B halo, whose satellites are somewhat closer to the 
main halo. 



6 



Table 3 

Reionization times of the main halos in the four models 

Model No. A B C D E F Box Mean 

1 9.2 9.9 12.0 9.3 8.5 8.1 9.31 

2 11.5 11.2 13.7 10.7 10.0 9.5 10.71 

3 12.6 12.4 16.3 12.0 11.8 10.9 11.97 

4 14.0 14.3 17.1 13.5 13.6 12.5 13.50 



galaxies. Our satellite properties are in good agreement 
with the observed values, and recover the observed rela- 
tion between M and ro.5. The points shown are for satel- 
lites with the Model 1 reionization history; driving reion- 
ization to earlier times yield the same relation, but fewer 
satellites overall, and fewer in the intermediate ro.5 (or 
luminosity) range. We conclude that our simple model 
is overall sufficient for reproducing the observed trends. 

4. RESULTS 

4.1. Spread in reionization redshifts 

A major novel factor in our work is that the reion- 
ization redshift of a halo is completely determined by 
its environment. This has two main consequences: first, 
that the six main halos will reionize at different times 
according to their locations in the Millennium-II box. 
But additionally, the subhalos of each of these main ha- 
los are not required to reionize at the same time as their 
main halo - their reionization times are also set by their 
particular path through the main box. Here, we explore 
the distribution in reionization redshifts for the subhalos, 
and how this infiuences the resulting number of potential 
satellites. 

Table[3]summarizes the first effect - that the six Aquar- 
ius host halos do not reionize at the same times in the 
box. In particular, the F halo consistently reionizes later 
than the box mean, while the C halo reionizes signif- 
icantly earlier. The other four halos are closer to the 
mean of the box (and thus slightly earlier than the me- 
dian; also see Table [2|). 

As for the second effect: Fig. |6] shows the distribution 
of satellite reionization redshifts for the six halos, here in 
the case of scenario 1. (The distributions for the other 
three models show similar characteristics, but with fewer 
satellites and narrowing distributions as the models be- 
come more efficient.) Several interesting trends are evi- 
dent: First, there is a distribution - while some of them 
are strongly peaked (C, D and F in particular), all show 
a range of redshifts. Second, the peak of the distribu- 
tion does not necessarily correspond to the reionization 
redshift of their host halo (the dotted lines) . 

The C halo is the most extreme in this regard - the 
peak of the satellite distribution is near the mean of the 
box, while the host reionized significantly earlier. Essen- 
tially, this happens because the location of the C halo is 
near the edge of one of these early-reionized regions in 
the large box (cyan in Fig. [1]), responsible for the early 
reionization of the C halo itself. Many of its satellites, 
however, fall in from a different direction than the bubble 
front, and end up with a reionization distribution peak- 
ing around the box mean. Note that Fig. [5] shows only 
systems that survive as self-bound satellites. If we were 
to also include systems that are shredded and are part 
of the diffuse main halo by redshift zero, we would see 



a larger number of subhalos reionizing around the same 
time as the C main halo. 

Having established that there is indeed a distribution 
of satellite reionization redshifts, the next question is 
whether this leads to a different satellite population be- 
tween patchy and instantaneous models. Fig. [7] consid- 
ers the number of satellites produced in each model for 
different sets of assumptions: either patchy or instanta- 
neous, and a suppression Vmax of either 30 or 50 km/s. 
We see that regardless of cutoffs and patchy vs. instan- 
taneous, as reionization is pushed to earlier times, there 
is a strong decline in the number of satellites. This over- 
all trend is similar to what was found in previous studie s 
(jMuiioz et al.|[200l lTumfinson|[2010l: iBusha et al.llMol) . 
For a given reionization model, we do see additional 
differences, however. For either choice of suppression 
Vmax, the instantaneous reionization models systemati- 
cally make 10-20% more satellites than the corresponding 
patchy models; in the more extended models this effect 
dominates the difference between stronger and weaker 
suppression. As the reionization sources become more 
efficient and reionization occurs earlier, the patchiness 
is less important, and the variation in suppression Vmax 
becomes the dominant effect in this scenario. 

Note that the numbers shown in Fig. [7] are averages; 
as shown in Fig. IH] there is a substantial spread between 
the actual number of satellites between the different ha- 
los. The error bar indicates the typical rms, ± 15 
satellites. However, halos A-E all show more satellite 
candidates in the instantaneous scenarios, while the F 
halo has about 10% fewer. For five out of the six ha- 
los, then, the main effect of patchy reionization is that 
the number of satellite candidates is reduced by 10-20% 
compared to the instantaneous models. This can be un- 
derstood in the context of Fig. [6l where the F halo is the 
only halo whose majority of satellites reionize later than 
the mean. 

In practice, the difference between the patchy and in- 
stantaneous case in Fig. [7] is due to a combination of 
two effects: that the global mean is not necessarily the 
best instantaneous redshift to choose for a given halo, as 
well as not capturing the spread in reionization times. 
However, we do not find that using e.g. the host halo 
reionization redshift, or the median satellite reionization 
redshift, leads to better agreement in each case. While 
we could recover the number of satellites in each case by 
tuning the instantaneous reionization redshift, the value 
of such an exercise is limited: the only way to determine 
that appropriate redshift would be to already know the 
result from the patchy model. As such, this illustrates 
the main problem with using instantaneous and/or spa- 
tially homogeneous reionization models. Depending on 
the application, the fact that there is a distribution of 
satellite reionization redshifts may be important to the 
results - and this distribution cannot be captured in a 
homogeneous model, nor can the appropriate instanta- 
neous approximation be recovered without knowing the 
patchy result. 

4.2. Luminosity Functions 

Having explored the reionization distributions in our 
simulations, we here discuss the resulting luminosity 
function of our simulated satellites. In order to com- 
pare them to the population of observed dwarfs, how- 
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Figure 5. Mass within the estimated half-Hght radius versus half-Hght radius for the simulated sateUites (red diamonds) for the six halos, 
in reionization scenario 1, with the higher suppression {vmax = 50 km s"'^). For the simulation points, the half-light radius is estimated 
from the luminosity via an empirical relation, and the mass is measured by fitting a NFW profile. Black triangles show the observed values 
in the Milky Way classical dSphs, while the crosses show the values for the UFDs. 
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Figure 6. Distributions of reionization redshifts of the satellite 
candidates of the six halos, here shown for Model 1 and the higher 
(vmax = 50 km s~^) suppression threshold. The dashed lines cor- 
respond to the mean reionization redshift across the entire box in 
Model 1, while the dotted lines show the reionization redshifts of 
the six main halos. Note the tails to higher redshifts. 



ever, we first need to correct the obser vational sam- 
ple of dwarf galaxies (data ta ken from iMartin et al.l 
l2008HWadepuhl fc Springe]|l2011[ ) for incompleteness and 
biases to establish a lumi nosity function . We fol- 
low the metho d outl ined in iKoposov et al.l ([2008) and 
iToUerud et al.l ()2008[ ). which employed models for the 
completeness of the SDSS Data Release 5 (DR5) to quan- 
tify the detection efficiency and thus adjust the observed 
luminosity function. In addition, we only include the 
simulated satellites that would be detectable according 
to the DR5 criteria: 

W = 10(-«-«^^-5-23)/3 Mpc, (2) 

where foRb = 0.194 is the fraction of the sky covered by 
DR5. That is, simulated satellites with r > rmaxiMy) 
are not included in the plots in this section. In each case, 
fmax is calculated with respect to a point 8 kpc from the 
center of the main halo. 

Figure [5] shows the resulting cumulative luminosity 
function for the simulated satellite population of all six 
halos, for the four patchy reionization models and with 
Vmax = 50 km s^^. The luminosities shown are those ob- 
tained with the abundance matching technique (eqn. [TJ , 
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Figure 7. Mean number of satellite candidates for the six halos, 
for our four different reionization models, with varying assump- 
tions; Either strong suppression at reionization, affecting all halos 
with Vmax < 50 km s~ , or a weaker suppression at 30 km s^^. In 
addition to the patchy cases, we also show calculations where reion- 
ization instead happened instantaneously at the mean redshift of 
each patchy model. While the trends with reionization model and 
suppression Vmax are similar for the six halos, there is substantial 
halo-halo scatter in the actual number of satellite candidates. This 
is indicated by the error bar (typical rms ~ 15 satellites). 



but the resulting luminosity functions are very similar 
to those using our simple SPS-based model. (However, 
there is substantial scatter between the methods on a 
subhalo-subhalo basis, i.e., rms ~ 1 mag.) The dotted 
line shows the resulting luminosity function when apply- 
ing the abundance matching prescription to all subhalos 
in each simulation that lie within the virial radius by 
redshift z — 0, without assuming any suppression due to 
baryonic processes. As such, it gives a good illustration 
of the missing satellite problem: without some process 
such as reionization to suppress star formation and re- 
duce the number of luminous satellites, the six Aquarius 
halos would predict between 500-1000 satellites around 
Milky Way-like galaxies. 

We see that while satellites brighter than My ~ — 10 
are insensitive to the reionization model, the total num- 
ber and the number of faint satellites change dramati- 
cally with more effi cient reionization. This is consistent 
with the results of iFont et aD (|2011l ) , who found that 
photoheating is mainly important for suppressing star 
formation in satellites analogous to ultra-faint dwarfs. In 
our simulations, only Model 1, and in some cases Model 
2, can match the faint end of the luminosity function. In 
fact, there is a factor of 3-4 difference between the to- 
tal number of satellites predicted between the different 
reionization models for a given halo. This behavior may 
ultimately allow for constraining the Milky Way reion- 
ization epoch with the faintest satellites. 

Another striking feature of Fig. |S] is the scatter be- 
tween the six halos - there is a factor of ^ 3 difference 
between the B and D halos, which have the fewest and 
most satellites, respectively. Thus, the halo-halo scat- 
ter is of the same order as the differences in the impact 
of the various reionization models on the faint end. In 
order to break this degeneracy, it will be necessary to un- 
derstand the general distribution of substructure around 
Milky Way-like halos, requiring a statistical sample of 



halos at sufficient resolution. Since the bright end of the 
luminosity function is not affected by the reionization his- 
tory, however, it may be possible to break the degeneracy 
by constraining the halo-halo scatter using the brighter 
satellites, and use the fainter satellites to constrain the 
reionization epoch of the Milky Way. 

We find that only the F halo hosts a satellite as bright 
as the Magellanic Clouds. This is consistent, however, 
with recent studies based on SDSS DR7 on the abun- 
dance of close, bright satellites around Milky Way-like 
gala xies (Liu ct al. 2011; Boylan-Kolchin ct al. 2011 a; 
iB^ha et al.ll2011l : iGuo et al.ir201ll: iToUerud et alllMTl . 
and the possibility that the Clouds have entered th e 
Milky Way halo only recently (jBesla et al.l 120071 12010D . 
Hence, we do not consider the mismatch at the bright- 
est magnitudes to be problematic, especially given the 
simplicity of our star formation model. 

Similar to the brighter satellites, the total stellar 
masses in the accreted halo of shredded, infalling satel- 
lites do not change significantly wit h reionization mode l. 
This is consistent with the results of lCooper et al.l ()201Clf) , 
who studied in detail the building of accreted stellar ha- 
los in the Aquarius simulations. They found that the sig- 
nificant progenitors are similar in mass to the brightest 
Milky Way dwarf spheroidals, with each halo assembled 
from less than five such objects. 

To highlight some of the differences between the vari- 
ous model assumptions. Fig. 12] shows the luminosity func- 
tion of the D halo satellites for some parameter choices. 
The upper panels show the results using patchy reion- 
ization and Vmax — 50 km s^^, instantaneous reioniza- 
tion and Vmax = 50 km s^^, and patchy reionization and 
I'maa; = 30 km s^^ respectively. The three lower pan- 
els show the same, but with luminosities determined by 
abundance matching (eqn. [T]). 

While the models with post-reionization star formation 
suppressed at Vmax = 50 km s~^ give a good match to 
the observed luminosity function, when we decrease the 
suppression threshold to 30 km s^^ none of the models fit 
the bias-corrected luminosity function well, despite the 
overall number of satellites only being slightly increased. 
This holds true both for the abundance matching-derived 
luminosities, and the best-ffi SFR/SPS model. In par- 
ticular, even if the star formation efficiency is decreased 
to give the correct number of satellites overall, the shape 
still does not match the observed luminosity function. It 
is possible, however, that this lower suppression model 
could still be a good match, if we ex plicitly consid er other 
processes like supernova feedback (|Li et al.ll20H . 

Compared to patchy reionization, the instantaneous 
models generally yield more satellites. While only the 
D halo is shown here as an illustration, as discussed in 
Section [4Tl this holds true for all of the Aquarius halos, 
except for the F halo. It is interesting to note, however, 
that the differences between different reionization mod- 
els, luminosity prescriptions, etc., are in general domi- 
nated by the halo-halo scatter. That is, the differences 
in the satellite populations of the six halos for any given 
model are as large as the differences between the range 
of models we are exploring. As an example. Figure [10] 
shows the satellite luminosity functions for all six halos 
overplotted, here for reionization model 1 and the SFR 
luminosity prescription. While the B halo, which is the 
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Figure 8. Cumulative luminosity functions of the simulated satellite populations of the six halos, using the abundance matching technique 
to assign luminosities. The colored curves show the results for the four patchy reionization models, with strong suppression at reionization 
[v-max = 50 km s""*^). The black points with error bars show the observed Milky Way satellite luminosity function, applying the completeness 
correction of Koposov et al. (2008) and ToUcrud ct al. (2008). The dotted line shows what the luminosity function for all the subhalos 
would be if we ignored all baryonic effects, and used the same abundance matching extrapolation to assign luminosities to all subhalos in 
each case. 



clear outlier, is also the lowest mass halo in this sim- 
ulation suite, the uncertainty of where the Milky Way 
itself fits in this distribution makes direct comparisons 
to constrain models difficult. 

5. DISCUSSION AND CONCLUSIONS 

In our models, we find that the population of faint 
satellites is highly sensitive to th e reionization mode l. 
This agree s with the results of iBusha et al.l (|2010l ): 
iFont et all (j2011[) . but is at odds with several earlier 
studi es in the literature ()Somervillell2002tlKravtsov et al.l 
|2004[) . The difference is largely due to two factors. First, 
the previous work that focused on the range of the clas- 
sical dwarf galaxies of the Milky Way found little reion- 
ization dependence; our results agree to the extent that 
the bright end of our luminosity function (classical dSph 
regime with 10^ < L < 10* L0) does not change either 
with large scale (different reionization model) or small 
scale (patchy versus instantaneous) changes in reioniza- 
tion. The difference at the faint end is explained by the 
various assumptions about how reionization influences 
the available supply of cold gas in small halos ~ our model 
assumes an abrupt end to star formation in halos smaller 



than a cutoff Vmax , whereas the filtering mass formalism 
will lead to a gradual suppression in star formation, and 
a much weaker reionization dependence. In that sense, 
we have chosen a maximum effect, and the results of our 
study can be viewed as an upper limit to the impact of 
patchiness of reionization. 

We also note that since all but the largest subha- 
los in our simulation finish their star formation around 
reionization, our results are consistent with some ultra- 
faint dwarf galaxies having formed around the epoch 
of reionization in ato mic cooling halos (but see e.g., 
IBovih fc Ricottill20Q9l ). Given the cutoff we set at the 
atomic cooling threshold, our model does not predict ob- 
jects much fainter than the faintest dwarf galaxies al- 
ready observed, although this could potentially change if 
II2 cooling at high redshift were included. 

Our main conclusions can then be summarized as fol- 
lows: 

• Patchy reionization models yield a distribution in 
satellite reionization redshifts, which leads to a 10- 
20% difference in the number of satellites they pre- 
dict compared to instantaneous models. This effect 
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Figure 9. Comparisons of ttie different models, fiere illustrated by the resulting luminosity function for satellites of the D halo. Upper 
left: same as Fig. |H] the patchy reionization model, with luminosities assigned by the SPS method, for Vmax = 50 km s"'^. Upper middle: 
instantaneous reionization at the mean reionization redshift of the model. Upper right: patchy, with threshold for suppressing star formation 
lowered to Vmax = 30 km . Lower panels: same as upper, but with luminosities assigned by the abundance matching technique, rather 
than the SPS method. 
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Figure 10. Illustrating halo-halo scatter: Overplotting the lumi- 
nosity functions of the satellite populations of all the six halos, here 
shown for reionization model 1 (yellow line in Figs. |8] and |9ll. 



can result in more or fewer satellites depending on 
the overall environment and merger history of the 
halo, and as such there is no general way to correct 



for it in a homogeneous model. 

• The overall number of satellites depends sensitively 
on the reionization model - we find a factor of 3-4 
difference between the earliest and latest reioniza- 
tion model. The difference is entirely in the number 
of fainter satellites, while the brighter end of the 
luminosity function corresponding to the classical 
dSph regime is unaffected. 

• There are large halo-halo variations in the number 
of luminous satellites for a Milky Way sized halo - 
we find a factor of 2-3 difference between the halo 
with fewest and most satellites, comparable to the 
spread between the different reionization models, 
lience, in order to use the Milky Way satellites to 
constrain reionization, a much larger sample of halo 
simulations is required. Given the spread between 
the six halos in the Aquarius sample, the direct 
comparison of one or a few halos to the observed 
Milky Way will only yield conclusive results once 
we have a more detailed understanding of how our 
Galaxy fits into this distribution. 

Progress can be made in several directions. On 
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the observational side, future and on going surveys like 
SkyMapper PanSTARRS and LSST (IKeller et alJl2007l: 
iKaiser et al.l[2002l: llvezic et al.|[2008l ) will cover the entire 
sky at the depth of SDSS or deeper, reducing the uncer- 
tainty in the number of observed dwarf galaxies due to 
the incomplete sky coverage of SDSS. Follow-up studies 
of new and existing dwarf galaxies can also provide more 
observational constraints, which can help discriminate 
between the different theoretical models that at present 
all fit the observed luminosity function. 

On the simulation side, the present work has shown 
that in a six-halo sample, there is a large degeneracy be- 
tween the significant halo-to-halo scatter and the results 
of effects like reionization. This level of cosmic variance 
will need to be quantified more accurately with a larger 
set of simulated halos. Only then can the origin and evo- 
lution of the faintest galaxies in the universe, as well as 
their role in building larger galaxies be more thoroughly 
understood. 
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